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1. Introduction 

Several types of multifractional Gaussian processes have been studied, including the 
processes introduced by Levy-Vehel and Peltier [12] and Bcnassi, Jaffard and Roux [6], 
known as the moving average and harmonizablc versions of multifractional Brownian mo- 
tion, and the processes introduced by Benassi, Cohen and Istas [2] and by Surgailis [16]. 
These processes are usually defined by replacing, in certain representations of fractional 
Brownian motion (fBm), the Hurst parameter H by a Hurst function h, that is, a real 
function of the time parameter with values in (0,1). Generalizations to stable processes 
have also been studied [10, 14, 15]. The stochastic properties and the regularity of the 
trajectories of these processes can be characterized in the same terms as for fBm, in 
particular, by the property of local asymptotic self-similarity or by the pointwise Holder 
exponent, which is constant and equal to H for fBm, but which varies in time following 
h in the multifractional case. 

The aim of this article is to study a Volterra-type multifractional Gaussian process 
Bh which fits into the framework of stochastic calculus, and its local time. In fact, the 
stochastic calculus, including the Ito formula and stochastic differential equations, is 
now well established for fBm and stochastic integrals have been defined in the Malliavin 
sense or by means of Wick products. The stochastic calculus for multifractional Gaus- 
sian processes has not yet been developed explicitly, but important elements, including 
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a divergence integral and an Ito formula, have been proven in [1] for a class of Gaus- 
sian processes admitting a kernel representation with respect to Brownian motion. The 
process Bh we study in this article is defined by allowing for a Hurst function in the 
kernel representation of fBm and it belongs to the class of Gaussian processes studied 
in [1]. The Ito formula in [1] is applied to obtain a representation of the local time for 
which, as is usually the case in Tanaka-likc formulas, the occupation measure is not the 
Lebesgue measure, but the quadratic variation process. However, contrary to fBm, the 
interpretation in our case is more delicate because the quadratic variation process is not 
necessarily increasing (but of bounded variation). We compare this local time to the local 
time with respect to the Lebesgue measure, which we obtain, together with the regularity 
properties of its trajectories, by the Fourier analytic approach initiated by Berman [3] 
based on the notion of local non-determinism (LND). 

It is well known that the cases where the Hurst parameter or the Hurst function is 
<l/2 (resp., >l/2) have to be treated separately (see [1] and [11]). We are interested in 
the case where the Hurst function h takes values in (1/2,1). In fact, the Volterra-type 
process we study here is defined only for this case, which is appropriate for long memory 
applications [13]. The article is organized as follows. In Section 2, stochastic properties 
and the regularity of the trajectories of Bh are proved. The form of the covariancc 
function (Proposition 2) shows, in particular, that Bh differs from the harmonizable 
and moving average multifractional Brownian motions. The continuity of the trajectories 
of Bh is obtained by classical criteria from estimates for the second order moments of 
the increments (Proposition 3). The lass property is proved in Proposition 5; it implies 
that the pointwise Holder exponent of Bh is equal to h (Proposition 6). In Section 3, 
it is shown that Bh satisfies the property labeled (K4) in [1] and that the quadratic 
variation of Bh is of bounded variation. Therefore, the divergence integral and the Ito 
formula developed in [1] (recalled briefly in this section) hold for Bh- Section 4 is devoted 
to the local time of Bh- Its existence and square intcgrability with respect to the space 
variable follow from a classical criterion due to Berman [3] (Proposition 11). A Tanaka- 
type formula with the divergence integral of Section 3 is then given (Theorem 12), where 
the occupation measure of the local time is the quadratic variation of Bh- Finally, the 
local time with respect to the Lebesgue measure and its regularity in the space and time 
variables are given in Theorem 16. The proofs are based on the LND property, which is 
shown to hold for Bh in Proposition 14. 

2. A Volterra-type multifractional Gaussian process 

It is well known that the fractional Brownian motion Bjj with (fixed) Hurst parameter 
H G (1/2, 1) can be represented for any t > as 
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where 



K H (t,u) =u 



= ,,V2-H 



(y-u) H -V\j H -^dy 



(1) 



and W(dy) is a Gaussian measure. Let a and b be two real numbers satisfying 1/2 < 
a<b < 1. Throughout the paper, we consider a function /i:M— >• [a, 6]. We assume that 
this function is j3- Holder with sup h < /?. Define the centered Gaussian process = 
{B h (t),t>0} by 



B h (t) = B h{t) (t) = / if h(t) (t,u)W(du) 



where 



K h {t)(t,u) =u 



_„l/2-fc(t) 



(2) 



Before establishing properties of -B^, we state a lemma regarding an estimate on Kh 
that we use throughout the paper. 

Lemma 1. For every T > 0, there exists a function G L ((0,T],R_)_) smc/i i/iai /or 
every s G (0, T], 



sup 

Ae[o,6],te(o,T] 



5 r, / \ 

—K\(t,s) 



< *t(s). 



Proof. We have 



+ s 1/2 ~ X f (y - S ) A ~ 3/ V- 1/2 (log(y - s) + logy)dy. 



Then, for every T > 0, there exists a constant C a bT such that for every A £ [a, 6] and 

se(o,T], 



<C Q , 6 , T (lV|l0g S |) S 1 /2-6 =: $ T ( s ) ) 



which concludes the proof. 

The following proposition gives the covariance of this process. 



□ 



Proposition 2. Let X = {X(t, A), t > 0, A G (1/2, 1)} &e i/ie two-parameter process given 
byX{t,X) = f*K\(t,u)W(du). Then 



i[X(t,X)X(s,X')} = dy dzp(y,z,X,X')\y-z\ 
Jo Jo 



x+x'-2( v 

Z 



A — A 
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where 

0(y, z, A, A') = 0(2 - A - A', A' - l/2)l {y>z} + 0(2 - A - A', A - l/2)l {y<z} 
and 0(a,b) (a,6>0) is the beta function. In particular, 



with 



E[X(t,X) 2 } = / dy / dz/3(2-2A,A-l/2)|y-z| 2A - 2 - 

JO JO C A 



2 = 2ttA(A - 1/2) 3 

A r(2 -2A)r(A + 1/2) 2 sin(7t(A- 1/2))' 1 ' 



Proof. We have 

E[X(t,X)X(s,X')] 



tAs 

K x (t,u)Ky(s,u)du 







tAs 



1XX ' {Iu iy ~ u)X ^' 2yX ^' 2 dy ) (Jj Z ~ «) A '~ 3/2zA '~ 1/2 dz ) <±u 

= dy dV" 1/2 ^ A '" 1/2 / w 1 - A - V (2/- U ) A - 3/2 (z-w) V - 3/2 du. 
■/o Jo Jo 

We fix y > z and calculate the following integral by making the successive substitutions 
u = vz, w = (y — vz) /(1 — v), t = w/y and s = 1/t: 



(y-u) x -^ 2 (z-u) x '-^ 2 du 

= z 1 ' 2 ^ f \'(y - vz) x -V 2 (\ - v) x '~ 3 ' 2 dv 
Jo 

= z 1 ^ x (y - z) x+x '- 2 / (w - y) 1 -^ w X ~ 3 / 2 dw 
Jy 

/ + oo 
(t-ly-^'^-^dt 

= z l ' 2 - x y l / 2 - x ' [y - z) x+x '- 2 0{2 - A - A', A' - 1/2). 



In the same way, for y < z, we obtain 



u 1 - x - x 







(y-u) x -V 2 (z-u) x '- 3 / 2 du 

z l/2-\ y l/»-\' (z _ y )A+A'- 2/?(2 _X-X',X- 1/2). 
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This concludes the proof. □ 
In the sequel, we need the estimates which we establish in the following proposition. 

Proposition 3. The process X satisfies the following estimates: 

(a) for all s and t>0 

E{(X(t,\)~X( S ,\)) 2 }=c- 2 \t-s\ 2X ; (4) 

(b) for every T > 0, there exists a constant Ct > such that for every t e [0,T] and 
every X and X' £ [a, b], 

E[(X(t, A) - X(t, A')) 2 ] < C T \X - A'| 2 . (5) 

Proof, (a) For every A £ [a, b] , the process X(-,X):t>-t X(t, A) is a fractional Brownian 
motion with variance c^" 2 , so we deduce (4). 
(b) We have 



E[{X(t,X)-X(t,X')Y}= / (K x {t,u)-K x >(t,u)ydu. 

Jo 

There exists £ = £(A, A') 6 [min{A, A'}, max{A, A'}] such that 
K x (t,u)-K x ,(t,u) = (X-X') 



d 

—K x {t,u) 



A=£=£(A,A') 



Then, thanks to Lemma 1, for every t, X and A', we obtain 

E{(X(t, A) - X(t, A')) 2 ] < |A - A'| 2 f ($ T (w)) 2 du, 

Jo 

which concludes the proof. □ 

From the last proposition, we can deduce the continuity of Bh- 

Corollary 4. The process Bh defined above has continuous trajectories. 

Proof. We deduce from Proposition 3 and the assumption sup h < j3 that for every s 
and t in a compact interval [0,T] such that \t — s\ < 1, we have 

E[(B h (t) - B h (s)) 2 } = E[(X(t, h(t)) - X(s, h(s))f] 
<2E[(X(t,h(t))-X(t,h(s))) 2 } 
+ 2E\(X(t,h{s))-X( S ,h(s))f} 
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<2C T \h(t)-h(s)\ 2 + 2clf s) \t-s\ 2h ^ 
<2C T \t~s\ 2f} + 2 sup (cf)\t-s\ 2a . 

\e[a,b] 

Then E[(Bh{t) — Bh(s)) 2 ]/\t — s\ 2a is bounded and since Bh is Gaussian, we can deduce 
from [5] its continuity. □ 

We now deal with the local self-similarity property of Bh- 

Proposition 5. The process Bh is locally self-similar. More precisely, for every t, we 
have the following convergence in distribution: 

( B h {t + eu)-B h (t) \ , 
l™o{ jHt) ) u>o = ( C h<t) B fc(*)(«))u>0> 

where lim £ _>o stands for the limit in distribution in the space of continuous functions 
endowed with the uniform norm on every compact set. 

Proof. Let us start by proving the convergence of the finite-dimensional distribution. 
Because Bh is Gaussian, it suffices to show the convergence of the second order moments. 
We then can write, for every u and v, 

~[{B h (t + eu) - Bh(t))(B h (t + ev) - B h (t))} 



£ 2h(t 

1 



{h(e)+I a (e)+h(e)+h(e)), 



£ 2h{t 

where 

h(e) = E[(X(t + eu,h{t)) - X(t,h(t)))(X{t + ev,h(t)) - X(t,h{t)))}, 

h{e) = E[(X(t + eu, h(t + eu)) - X(t + eu, h(t)))(X(t + ev, h(t)) - X(t, h(t)))}, 

h(e) = E[(X(t + eu, h{t)) - X(t, h(t)))(X(t + ev, h(t + ev)) - X(t + ev, h(t)))], 

h(e) = E[(X(t + eu, h(t + eu)) - X(t + eu, h(t)))(X(t + ev, h(t + ev)) - X(t + ev, h(t)))}. 

By the self-similarity of the fractional Brownian motion and the stationarity of its incre- 
ments, we obtain 

^M*) = ^(\-\ 2h{t) + M 2hit) - 1" - y\ 2h(t) )- (6) 

£ h(t) 

Then, because of the Cauchy-Schwarz inequality, it is enough to prove 

lim -5^7jrE[(A'(t + eu, h(t + eu)) - X(t + eu, h(t))f] = (7) 

e— >0 £*> l V') 
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to get 



lira 



1 



e-X) £ 2 M* 

Thus, using Lemma 1, we obtain 
1 



(h(e) + h(e) + h(e))=0. 



(8) 



?2h(t) 



E[{X{t + eu, h(t + eu)) - X(t + eu, h(t))Y 



1 



-2h{t) 



t-\-eu 



(Kh(t+eu)(t + eu,s)-K h{t - ) (t + eu,s)) ds 



< 



(h(t + eu) - h{i)f 

F 2h(t) 



2 rt+eu 



sup 

H 



9 „ , 
—K H {t + eu,s) 



ft+u 

Jo 



It now remains to prove tightness in the space of continuous functions endowed with the 
uniform norm. We also consider T > such that t, t + eu and t + ev £ [0,T] for all e. 
Performing similar calculations as in the proof of Corollary 4, we get that there exist 
Ct > such that 



E 



B h (t + eu)-B h (t) B h (t + sv)-B h (t) 



E[(B h (t + eu) - B h (t + ev)) 2 } 



< 



z 2h(t] 

Ct 



\eu — ev 



2h{t+ev) 



£ 2h(t) ■ 

= C T E 2{h{t+ev) ~ h{t)) \u - V \ 2h ^+ £V \ 

Since h is /3-H61der, it follows that £ 2 ( h (t+^)-h(t)) - 1S um f orm iy bounded. Moreover, h(t- 
ev) > a, thus 



B h (t + eu)-B h (t) B h {t + sv)-B h {t) 



£ h(t) 

This completes the proof. 



Mt) 



2t 



< Ct \u — v I 



□ 



It is classical to deduce pointwise Holder continuity from local self-similarity [2] . We re- 
call that the pointwise Holder continuity of a function / is characterized by the pointwise 
Holder exponent aj(to) defined at each point to as 
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Proposition 6. For every to GR + , the pointwise Holder exponent ctB h (to) of Bh is 
almost surely equal to h(to). 

Proof. We deduce from Proposition 5 and [2] that ctB h (to) < h(to)- We now prove that 
olb h (*o) > h(to). By Proposition 3, there exists C > such that for every e > and every 
s and t € [to — e, to + e] such that \t — s\ < 1, wc have 

E[(B h {t) - B h (s)f] <C\t- ,s| 2inf [*o--'o+ei \ 

By using the fact that -B^ is Gaussian and applying Kolmogorov's theorem [5], we get 
that lim t _>.t \B h (t) - B h (t )\/\t - t Q \ a = for every a < inf [to _ £jto+£ ] h. This holds for 
every e > 0, so by continuity of h, we have lim t ^ to \Bh(t) — Bh(to)\/\t — to\ a = for 
every a < h(to). We can deduce that aB h (to) > h(to) and hence a.B h {to) = h(to). □ 

3. Stochastic calculus for Bh 

The aim of this section is to apply the stochastic calculus developed by Alos, Mazet and 
Nualart in [1] to get a stochastic integral for Bh and an Ito formula. We recall that in 
[1], the following hypothesis, called (K4), appears for regular kernels: 

• Hypothesis (K4). For all s € [0,T), K(-,s) has bounded variation on the interval 
(s,T] and J Q T \K\{{s 7 T],s) 2 ds < oo (where, for all sE [0,T), \K\((s,T],s) denotes 
the total variation of K(-,s) on (s,T]). 

Lemma 7. Suppose that h is of bounded variation on (s,T] for all s € [0,T). XTien (K4) 
/io/rfs /or (t, s) h-> Khn^t, s) defined by (2). 

Proof. Let 

n 

Var? T1 (-,s) = sup V|^ (ti) (^, s) - K h{t . _ l) (t i _i,s)| 

to=8<ii<-<t„=T^ 

and suppose, without restriction of generality, that T = 1. Then 

|^(t i) (ti,s)-if fe( t i _ l) (t i _ 1 ,s)|</i(i) + 7 2 (i), 

where 

■MO = l^(ii)(*i' S ) — -K'ft.Cti) i , ») I 

and 

h{i) = \K h{u) (t t _ ll s) - K h ( t ._j(ti-i,s)\. 
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We have 



7 1 (i) = S V2-M^) / (y- S )h(U)-3/2 y H ti )-l/2 dy 



< s 



1/2-b 



Therefore, 



^2\ K h(t,)(ti,s) - K h ( ti) (ti-i,s)\ 



i=l 



<C(a)s l ' 2 - b . 

Regarding I 2 , we have 

n 

d 



(9) 



<^2\h(ti) - hiti-Jl sup 



<t<l,a<\<b 



dx 



K x (t,s) 



< Var h sup 

[0,T] s <t<l,a<\<b 



d_ 

dX 



K x (t,s) 



The proof of Lemma 1 implies that Var™ s T j (•, s) < C(l V | logsDs 1 / 2 fc , where the constant 
C > depends on h (but not on n). This implies that (K4) holds for Kh- □ 

Remark 8. Since /i is also supposed /3-H61der continuous for some f3 < 1, we will here- 
after suppose that h is Lipschitz continuous. Note that limtx^ K^-jit, u) = for all 
u>0. 

For simplicity, we write K instead of Kh, but keep in mind that (i, s) — > K(t, s) means 
{t,s) — > Kh{t){ti s) and that differentials with respect to t also act via h. 

In the sequel, we need the following proposition regarding the variance of B/, . 

Proposition 9. The variance s i— > R s =E[_B/,.(s) 2 ] is of bounded variation on (0,T]. 

Proof. We have 



K(s t+1 ,r) 2 dr- / if(s l; r) 2 dr 
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i=i 



<E / K{s l+ll rfdr + Y. 



i=i 



[K(s i+1 ,r) 2 - K( Sl ,r) 2 ]di 



The functions |J"T(s,r)lm s i(r)| are bounded by the square-integrable function fc(r) = 

\K(T,r)\ + \K |((r, T\,r). Hence, the first term above is upper bounded by k(r) 2 dr. 
The second term is upper bounded by 

E / \K\ ((«i, s i+1 ], r) \K(s i+1 ,r)+K( Si ,r)\ dr 



<2 / k{r)\K\((r,T],r)dr < oo. 
This concludes the proof. 



□ 



For any / e L 2 ([0,T]), let Kf be defined by (Kf)(t) = /* K(t,s)f(s)ds. Let £ be the 
set of step functions on [0,T] and let the operator K* be defined on 5 by (K*ip)(s) = 
ip(t)K(dt,s). Then K* is the adjoint of K. In fact, for a* g M, = Si < s 2 < • • ■ < 
s n+ i=T, y = J2i=i a i 1 (s i ,s i+1 ](s) and /ei 2 ([0,T]), we write 

n n 

= E ( s ) E «) - *)) 



i=l 



and 



/ ( j PTV)(s)/(s)d S = E^ / El(. <) . <+1 ](*)(«'(*j+i,*)--K'(»J,*))/(*)d« 
Jo j=1 Jo . =1 

= f>i[(#/)(*i+i)-(tf/)(a^ = f V {t){Kf)(dt). 

3=1 J ° 



(10) 



As usual, the reproducing kernel Hilbcrt space (RKHS) % is defined as the closure of the 
linear span of the indicator functions {l[o.t],i € [Oi^ll with respect to the scalar product 
(l[o,t],l[o,«])ft =E[B /l(t) (t)B /l(s) (s)] = R{t,s). By replacing /(s)ds in (10) by W(ds), we 
have, by (1), 



B h {^)= cp(t)B h (dt)= (K*^)(s)W(ds). 



Therefore, 



||^||^= E[i? /1 M 2 ] = ||KV||| 



2 ([0,T]) 



< 



I^WII^Kdt.a) 



d*=:|M&. 
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Let us denote by Hk the completion of £ with respect to the || • ||^-norm. Then Hk 
is continuously embedded in %. 

In order to define the stochastic integral with respect to Bh , let us denote by S the set 
of smooth cylindrical random variables of the form F = f(Bh(<pi), Bh{f2), ■ ■ ■ ,-B(,(^ n )), 
where n > 1, / G C£°(R") (/ and all its derivatives are bounded) and ipi, <p2, ■ ■ ■ , <p n 
Let us also denote by n^ 2 {U K ) the closure of {F G S : F G L 2 (£l,U K ),DF G L 2 {Q x 
H K ,H K )}- Then D 1 ' 2 ^) is included in the domain Dom(,5 Sh ) of the diver gence oper- 
ator of Bh, and the integral of u G Dom(8 Bh ) with respect to Bh is given by 



T r T pT 



5 Bh {u) = I u8B h = / (K*u)SW 



T 



u(r)K(dr, s) 



5W(s), 



where the last two integrals are the divergence integrals with respect to Brownian motion. 
Let us recall that the integral 5 Bh (u) is defined, for any u <G L 2 (n,"H), as the unique 
element in L 2 (Vi) which satisfies the duality relationship E(S Bh (u)F) ~E(DF, u)u for 
all F G S. 

The following Ito formula, due to Alos, Mazet and Nualart [1], will be applied in 
the next section. Let F be a function belonging to the class C 2 (R) and satisfying the 
condition 

max{\F(x)\, \F'(x)\, \F"(x)\} < ce xlxl \ 

where c and A are positive constants such that A < i(sup < t<T i?t) _1 . This implies that 

the process F'(B h ) belongs to B 1 - 2 ^). The integral f* F'{B h {s)) 5B h {s) is therefore 
well defined for all t G [0,T] and the following Ito-type formula holds ([1], Theorem 2): 

F{B h (t))=F(0)+ [ F'(B h (s))6B h (s) + l [ F"{B h (s))dR(s). (11) 
Jo 1 Jo 



4. Local time and Tanaka formula for Bh 

First, we prove, by means of a criterion for Gaussian processes due to Berman, that Bh 
has a local time with respect to the Lebesgue measure. We then derive a Tanaka-type 
formula from the Ito formula of Section 3 and show that Bh satisfies the LND property 
This implies continuity and Holder regularities in space and in time of the trajectories 
of the local time. 

Definition 1 0. For any Borel set C C R+ , the occupation measure mc of Bh on C is 
defined, for all Borel sets Act, by mc{A) = \{t G C, B h ^{t) G A}, where A is the 
Lebesgue measure on R + . If mc is absolutely continuous with respect to the Lebesgue 
measure on R ; the local time (or occupation density) of Bh on C is defined as the Radon- 
Nikodym derivative of mc and will be denoted by {L(C, x), x G R}. Sometimes we write 
L(t,x) instead of L([0,t],x). 
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This definition implies that the local time of Bh satisfies the occupation density formula 



g(t,B h (t))dt = / g{t,x)L{dt,x)dx (12) 

C Jcxl 

for all continuous functions with compact support j:(7xt->R. If g docs not depend 
explicitly on t, then we get the more classical occupation density formula where the 
right-hand side of (12) is replaced by J M g(x)L(C, x) dx. 

Proposition 11. The local time of exists P-a.s. on any interval [0,T] and is a 
square-integrable function of x. 

Proof. By [3], for any continuous and zero-mean Gaussian process {X t , t 6 [0,T]} with 
bounded covariancc function, the condition 



dsdt 



oJo y/E[\X t -X s \< 



<+oo (13) 



is sufficient for the local time of X to exist and to be a square-integrable function of x. 
For (s,t) away from the diagonal, we write, for s <t, 

E[\B h (t) - B h (s){ 2 } = J (K h{t) (t, u) - K h{s) {s, u)f du + J K h(t) (t, u) 2 du 

and deduce from (1) that the second term stays strictly positive as (s,t) varies in [0,i — 
e] x [0,T]. If (s,t) is close to the diagonal, say < t — s < e, then, by Proposition 5, 

E[\B h (t)-B h (s)\ 2 }~c-* t) (t-s) 2h W 

as s /~t and a direct calculation shows that (13) is satisfied. □ 

Let us now derive a Tanaka- type formula for Bh- Since the last term in the Ito formula 
(11) is an integral with respect to the variance function R and since R is not, in general, 
increasing, but only of finite variation, the trajectorial representation of the local time 
is more delicate than for Brownian motion or for fractional Brownian motion. In fact, 
on the time intervals where R is (strictly) increasing or decreasing, this formula gives an 
occupation density L related to L. 

Theorem 12. Suppose that h is continuously differentiable. Then, for all o£R, 

\B h (t)-a\ - \B h (s)-a\= J sign(B h (u) - a) 8B h (u) + L([s,t], a) 

P-a.s. , where L([s,t],a) = J"* i?'(u)L(du, a) and L is the local time with respect to the 
Lebesgue measure of Bh- On the time intervals [s,t] where R is strictly increasing (resp., 



1306 B. Boufoussi, M. Dozzi and R. Marty 

strictly decreasing), L (resp., —L) is the (positive) occupation density of Bh with respect 
to the measure induced by R. 

Remark 13. (a) No information on the local time of Bh can be deduced from the above 
formula if R'(u) = on the interval [s,i\. In fact, 

R'(u) = ±nB 2 h{u) (u)] = ^(^ (u) ) = 2(fc'(«)logu+ ^t^M 

and R'(u) = if h(u) = 1/logu £ (1/2, 1) on an interval (s,t). In this case, L([s,t],a) = 
for all a. 

(b) Tanaka formulas for fractional Brownian motion have been shown by several au- 
thors. We refer to the survey [8] for references. 

Proof of Theorem 12. For e > 0, let p e (x) = (2ne)~ 1/2 exp(-x 2 /(2e)). We apply the 
Ito formula of the previous section to the function 

f b (x)= f F'M&y, 



where 

F' e {x)=2l p s {y)dy-l. 



F' £ {x) then converges to sign(x) and F e (x) converges to \x\ as e — s- 0. By (11), for e > 
fixed, 



F e (B h {t) - a) = F e (-a) + / F' e (B h {r) - a)K(dr, s) SW S 

J ° ^ d4) 
Pe (B h (s)-a)dR(s). 

Note that, by (K4), the process {F' £ (B h (r) - a),r e [0,t]} belongs to L 2 (n,n K ) and 
therefore to Dom(S Bh ) . Or, equivalently, the process {/ s * F' e (Bh (r) — a)K~h (dr, s), s € [0, i] } 
belongs to Dom((5 lv ). Clearly, F £ (Bh(t) — a) converges to \Bh(i) — a\ in L 2 (S1) and F e (—a) 
converges to \a\ if e— > 0. Moreover, the process {F' e (Bh(r) — a),r £ [0,t]} converges, as 
e^0, to {sign(B h (r)-a),r€ [0,i]} in L 2 (n,H K )- In fact, by (K4), 

E f(f s \Fi(B h (r)-a)-siga(B h (r)-a)\\K h \(dr,s)\ ds 

is upper bounded independently of e and we may apply the dominated convergence 
theorem ([9], Lemma 1). Therefore, the last term in (14) converges in L 2 (£l). Let us 
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denote the limit by A". Therefore, for any continuous function / with compact support 
in R, 



Pe (B h (s)-a)dR(s) /(a) do 



(15) 



converges in L 1 (f2) to J A"/(a)da. In fact, the dominated convergence theorem applies 
because 



Jo 



E[p £ (B h (s) - a)} dR(s) 



f PR(s)+e(a)dR(s)< f s- b R'{s)ds<oo. 
Jo Jo 



But (15) also converges to L f(Bh(s)) dR(s) = L f(Bh{s))R' (s) ds, where we use the 
fact that R is diffcrentiablc if h is differcntiablc. Hence, 



f(B h (s))R'(s)ds = / A?/(a)da 



By the occupation density formula (12) applied to g(s,x) — f(x)R'(s), we get A" = 
J Q R'(s)L(ds, a) = L(t, a) for A-a.e. a, P-a.s. We can extend this to all a £ R since there 
exists a jointly continuous version of A, and therefore of L, as will be shown next by 
Berman's methods, which are independent of the Tanaka formula. □ 

We now state regularity properties in time and space of the trajectories of L. The 
regularity properties of L follow easily since L(t,x) — J Q R'(s)L(ds,x). In order to show 
the existence of a jointly continuous version of A, we recall the hypotheses introduced in 
[7] and show that they are satisfied for Bh- We recall them in terms of any real- valued 
separable random process {X(t),t £ [0,T]} with Borcl sample functions. 

• Hypothesis (A). There exist numbers po > and H £ (0, 1) and a positive function 
ip £ i x (R) such that for all A £ R and t, s £ [0,T], < \t — s\ < p , we have 



E 



X(t)-X(s) 

exp iA 



\t-s 



H 



<V(A). 



Hypothesis (A m ). There exist positive constants 8 and c (both eventually depend- 
ing on m > 2) such that for all ii,t2, • • • ,t m with =: to < t\ < ■ • ■ < t m < T and 
\t m — tx \ < S, we have 



aP^E«i(^(tj)-^M)) 



<niE[cxp(ic U ,(A(t 3 )-A(t J _ 1 )))]| 
i=i 



for all Ui,U2, ■ 



i £ ■ 



Hypothesis (A) is evidently satisfied for self-similar processes with stationary incre- 
ments: here, ip{\) = \K[e lXX ^]\. Hypothesis (A) is also closely related to asymptotic 
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self-similarity and, in fact, it holds for a fairly large class of processes (see Proposition 
15). 

If X has independent increments, (A m ) is trivially true for all m > 2. When the left- 
and the right-hand sides of the inequality are applied to the characteristic function of 
a Gaussian process X, we get the condition which is known in the literature as local 
non- determinism (LND). Loosely speaking, LND says how much dependence is allowed 
for the increments of the process if the local time has certain regularity properties. As a 
general rule, the trajectories of local time become more regular if the trajectories of the 
process become less regular. 

Proposition 14. For all T > 0, the process {Bh(t),t £ [0,T]} satisfies hypotheses (A) 
and (A m ) for all m>2 with H = maxo<t<T h(t). 

Proof. Let us start by showing (A) . By similar calculations as in the proof of Proposition 
5, for every [0,T], 

[ e -*W>(B h (t + e)-B h (t)) 2 ]--2- <e^ s ^ T] h f T+ \$ T+1 ( s )f ds . 

c h{t) Jo 

We can now conclude by applying Proposition 15. 

We now prove (A m ) for all m > 2. Wc show that Bh satisfies the LND property as 
introduced by Berman for Gaussian processes. For simplicity, we write B instead of Bh 
in this proof. Let t\ < t% < ■ • • < t m and let V m be the relative conditioning error given 

by 

_ Vax[B(l^)-B(t m - 1 )\B(t l ),...,B{t m -i)] 

V m — 



Wax[B{t m ) - B{t 
The Gaussian process is said to be LND if 

liminf V m > 0. (16) 

c\0+ 
0<t m -t t <c 

This condition means that a small increment of the process is not completely predictable 
on the basis of a finite number of observations from the immediate past. For Gaussian 
processes, Berman [4] has proven that (16) implies (A m ). More precisely, he has shown 
that if X satisfies (16) for all m > 2, then there exist constants C m > and S m > such 
that, for all u\,U2, ■ ■ ■ , u m £ R, 



Var 



Y^^iB^) ~ B(t 3 -i)] 
U=i 



3=1 



where to = and t% < • • • < t m are different and lie in an interval of length at most S m . 
This implies (A m ). In order to prove that B verifies (16), fix m > 2 and let < t\ < t^ < 
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< t m < T. By Proposition 3, 

VB,T[B(t m )-B(t m - 1 )}<C HT) (t m -t m ^ 1 ) 2h ^\ 



(17) 



where Ch(r) is a constant depending on T. We now study the conditional variance. 
Because of the definition of Bh in terms of W, we have the inclusion of cr-algebra 



We then get 



a{B{h), . . . , B(* m _i)} C a{W(t),te [0,i m _i]}. 

Var[fl(i m ) - B^-OlBfa), . . . , B(t m _i)] 
> Var[B(t m ) - B{t m - X )\W{t),t e [0,i m _ x ]] 
= Var[fl(t m )|W(t),te[0,t m _i]] 



(18) 



= Var 



#h(t TO )(*m,«)W( du ) 



(% m )(*m-' i )) d?X. 



We simplify the notation by letting t := t TO , s := i m _i and i? := h(t) and considering 
f*(KH(t 7 u)) 2 du. We recall that we have the relation 

/ (K„(t,u)) 2 du = E[(B H (t) - B H {s)f] - f (K H (t,u)-K H (s,u)) 2 du (19) 

Js JO 

so we study f Q (Kff(t, u) — Kh(s 7 u)) 2 du. We make the same substitutions as in the proof 
of Proposition 2 and thereby obtain 

{K H (t, u) - K H (s, u)f du 



.1-2H 



(y-u) H - 3 / 2 y H - 1 / 2 dy) du 



= dy dzy^^z"- 1 ' 2 / u^iy-u^-^iz-uf-^du 

Js Js Jo 



dy dz \y - z 



2H-2 



10 



(v-l) 1 ~ 2H v H -V 2 dv, 



where 



_ , . yz — zs _ 



yz - ys 



y>z "T l?/<2- 

yz — ys yz — zs 
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Then, using (19) and making the substitutions (y, z) — > (y + s, z + s) and (y, z) — > ((t 
s)y,(t-s)z), we get 



(K H (t,u)) 2 du = (t-s) 2H dy dz\y-z\ 2H - 2 

Jo Jo 



(v-l) 1 - 2H v H ^l 2 dv, 



where 



_ j . (t — s)yz + ys (t — s)yz + zs 

V 2 {s,t,y,z) = ) ( y y _ l y>z + -g- — l y <,. 



(t — s)yz + zs (t — s)yz + ys 



We can therefore deduce that 



where 



with 



liminf (t-s)' 211 [ (K H (t, u)) 2 du > X, 

(t-a)-*) J s 

1= dy dz\y-z\ 2b - 2 / (v^l) 1 - 2b v a - 3 ^ 2 dv>0 
Jo Jo J*3(y,z) 



(20) 



*3(j/)2) =SUp 

Combining (17), (18) and (20), we then get 



liminf V m >— 

0<t m -ti<c 



>o, 



which concludes the proof. 



(21) 



□ 



The beginning of the proof of Proposition 14 shows that the hypothesis (A) holds for a 
much larger class of processes. In fact, we have only used the assumption of the following 
proposition. 

Proposition 15. Let {X t ,t £ [0,T]} be a centered Gaussian process. Suppose that for 
some positive continuous functions f : [0, T] — > (0, 1) and g : [0, T] — > (0, oo), 



E[e- 2f{t) (X(t + e)-X(t)Y 



(22) 



uniformly in t. Hypothesis (A) then holds with H = sup/. 
Proof. Let us fix H = sup/. Because X is Gaussian and centered, we have 
X(t + e)-X(t) 



E 



cxp iA 



(23) 
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Because of (22), there exists £o such that for every e satisfying |e| < £o and for every t, 
we have 



E 



X(t + e)-X(t) 
7TW) 



c 



where c = inf 



[o,T] 9 



Besides, £ 2 /( t )- 2ff > l and thus 



(24) 



E 



X(t + e)-X(t) 



. £ 2f{t)-2H E 



X(t + e)-X(t) 



c 

> -. 

- 2 



(25) 



Combining (23) and (25), we then get, for every A, and t and s satisfying \t — s\ < eq, 
that 



E 



exp iA 



X(t + e)-X(t) 



tH 



< 



exp 



We then choose tp(X) = exp(— A 2 c/4) to conclude the proof. 



□ 



Let us now state a regularity result for the trajectories of L. The following theorem 
states that L(t,x) is Holder continuous in t of order 1 — H and Holder continuous in 
x of order , where H is the constant appearing in (A). For the proofs, we refer to 
[7], where it is shown that these regularities hold for any process starting at zero and 
satisfying (A) and (A m ) for all to > 2. Here, the function ip in (A) satisfies 



| u |C 1 - fl )/ ff ^( u )du<oo. 



| 14 | > 1 



Theorem 16. {Bh(t),t £ [0,T]} has a jointly continuous local time {L(t,x),(t,x) S 
[0,T] x K}. Moreover, for any compact set K C WL and any interval I C [0,T] with length 
less than po (the constant appearing in (A)): 

(i) ifQ<i<{l-H)/2H, then \L(I,x) - L(I,y)\ <r)\x - y\* for all x,yeK; 

(ii) if0<8< 1-H, then sup x g^- L(I , x) < , where r\ is a random variable, almost 
surely positive and finite, and \I\ is the length of I . 
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